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Abstract 

In lakes, trophic change and climate change shift the relationship between nutrients and physical factors, like temperature 
and photoperiod, and interactions between these factors should affect the growth of phytoplankton species differently. We 
therefore determined the relationship between P-limited specific growth rates and P-quota (biovolume basis) of 
Stephanodiscus minutulus and Nitzschia acicularis (diatoms) at or near light saturation in axenic, semi-continuous culture at 
10, 15 and 20 °C and at 6, 9 and 12 h d~^ photoperiod. Photoperiod treatments were performed at constant daily light 
exposure to allow comparison. Under these conditions, we also performed competition experiments and estimated relative 
P-uptake rates of the species. Temperature strongly affected P-limited growth rates and relative P uptake rates, whereas 
photoperiod only affected maximum growth rates. S. minutulus used internal P more efficiently than N. acicularis. N. 
acicularis was the superior competitor for P due to a higher relative uptake rate and its superiority increased with increasing 
temperature and photoperiod. S. minutulus conformed to the Droop relationship but N. acicularis did not. A model with a 
temperature-dependent normalised half-saturation coefficient adequately described the factor interactions of both species. 
The temperature dependence of the quota model reflected each species' specific adaptation to its ecological niche. The 
results demonstrate that increases in temperature or photoperiod can partially compensate for a decrease in P-quota under 
moderately limiting conditions, like during spring in temperate lakes. Thus warming may counteract de-eutrophication to 
some degree and a relative shift in growth factors can influence the phytoplankton species composition. 
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Introduction 

Interactions between nutrients and physical factors like 
temperature and light are important for phytoplankton growth. 
Nutrients and light are interdependent in surface waters because 
low nutrient levels restrict the phytoplankton biomass, and 
therefore typically coincide with clearer water. In turn, clearer 
water leads to higher underwater irradiances [1] and a longer 
efFective photoperiod due to a deeper euphotic depth [2]. 
Turbidity on the other hand increases absorbed radiation and 
therefore influences the water temperature and thermal structure 
[3,4], and water temperature affects nutrient cycling processes [5]. 
Therefore, in addition to seasonal cycles, the relationships between 
nutrients and physical factors can change, for example, due to 
global warming or eutrophication [5,6], with consequential efTects 
on phytoplankton communities that can be difficult to predict. For 
example, the dominance of filamentous cyanobacteria during 
spring in a shallow lake depended on the combined effects of 
winter temperature and the P:Si ratio [2] and the timing of the 
phytoplankton bloom was synergistically affected by water 
temperature and phosphorus supply [6]. The occurrence of 
several Aulacoseira (formerly Melosird) species could be differentiated 
along an irradiance/photoperiod - phosphorus gradient, suggest- 



ing trade-offs between these factors [7]. Understanding these 
effects may depend on how well we understand the physiological 
response of individual species. 

While Liebig's Law of the Minimum is assumed to apply to 
interactions between nutrients, this is not the case between 
physical factors like temperature, light and the photoperiod [8], or 
between nutrients and physical factors [9]. Temperature and 
photoperiod affect maximum specific growth rates of phytoplank- 
ton [10] and were shown to have species-specific interactive effects 
on nutrient-replete growth, where for example the response to 
changes in the photoperiod under constant daily fight exposure 
depended on temperature [11]. However, these interactions may 
be different under nutrient limitation because temperature and 
light not only affect growth rates, but also N- and P-quotas [12— 
14]. At the same time, light and temperature also influence 
nutrient uptake rates in a nutrient- and species-specific manner 
[15], whereby the temperature dependence of uptake is typically 
different to that of growth [16]. 

While the photoperiod was shown to have an important effect 
on spring phytoplankton [11], its interaction with phosphorus 
limitation is relatively unexplored. If phosphorus uptake rates 
differ in the dark and light [17], then there may be an interaction 
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between phosphorus and the photoperiod [18], particularly given 
the influence of light on phytoplankton stoichiometry [19]. 
Litchman et al. [20] showed that the combined effects of 
photoperiod and P-limitation were species-specific and greater 
than the sum of individual effects for several phytoplankton 
species, whereas Riegman and Mur [21] found a more either/or 
type of limitation between photoperiod and phosphorus for 
Pkmktothrix (formerly Oscillatorid) agardhii. Overall, it seems difficult 
to form a conclusive picture from these studies on photoperiod 
interactions with phosphorus, also because comparison is difficult 
when the irradiance and not the daily light exposure is held 
constant in daylength treatments. Under nutrient-replete condi- 
tions, light-saturated growth rates do not increase linearly with the 
photoperiod but show saturation characteristics [8,11,22-25]. 
Because of this nonlinear response, the effect of the photoperiod 
on growth rates can only be compared if the daily light supply is 
kept constant. 

Interactions between phosphorus and physical factors have 
become particularly relevant because recent research has stressed 
the importance of developing a mechanistic or biochemical basis 
for Droop's [26] quota model, which successfully describes 
phytoplankton growth rates, stoichiometry and the dynamics of 
limiting and non-limiting nutrients [27,28]. Droop's model 
[/^ = /i'm(l-Qo/Q)] relates the growth rate (jJ) to the nutrient 
quota (Q) in terms of the minimum quota (Qp) and the theoretical 
maximum growth rate at infinite quota Qo is the quota at 

which growth is zero and represents the amount of nutrient 
required for cell structure and machinery' [28]. Q,,, is the 
maximum quota at the real maximum growth rate when the 
nutrient is not limiting. The quota can exceed under luxury 
consumption, but this excess P is essentially stored as polypho- 
sphate and does not allow an immediate further increase of the 
growth rate [29]. The form of Droop's quota curve is fixed by the 
ratio Qc:Qm, "r expressed differently, such that the curve is half- 
saturated when the quota is double Qp. Qp seems to decrease with 
increasing temperature [13,30-32], although there are exceptions 
[33,34]. Furthermore, not only the maximum growth rate and Q, 
but also Qp:Qni [30] are all temperature dependent. Because the 
upper part of the quota curve is probably more important in 
deducing competitive advantage between species than the lower 
part near Qp [35], the Droop model may not be sufficient in 
characterising temperature or photoperiod interactions with P- 
limited growth. Of equal ff not greater importance than the 
growth-quota relationship are the nutrient uptake kinetics, and the 
feedback between uptake and quota [27]. The growth-quota 
relationship thus needs to be considered in conjunction with 
uptake. 

In this study we therefore investigate the effects of temperature 
and photoperiod on phosphorus-limited growth by comparing the 
parameters of the P-quota curve fitted to P-limited growth rates of 
two diatoms measured in semi-continuous culture under different 
temperatures and photoperiods. Due to the typical transition from 
limitation by physical factors to nutrient limitation during the 
course of spring, we assume interactions may play an important 
role here and accordingly, we chose two spring blooming diatoms 
as test species. We also examine the shape of the quota curve and 
the applicability of Droop's relationship. To enable comparison 
between photoperiod treatments, we adjusted the irradiance in 
experiments to compensate for different photoperiods and 
maintain a constant daily light supply [11,36-39]. Predicting the 
outcome of competition under nutrient limitation depends largely 
on obtaining accurate uptake measurements, which can be difficult 
[40,41]. As an alternative, we performed competition experiments 
to deduce the relative uptake affinity of the two species based on 



rates of competitive exclusion and the measured growth-quota 
kinetics. According to Healey [42], we assume that uptake affinity 
and competitive ability are proportional to the initial slope of the 
Michaelis-Menten uptake curve. Finally we extend an existing 
model that describes temperature-photoperiod-light interactions 
under nutrient replete t:onditions [1 1] to include phosphorus 
limitation. Our hypothesis is that both temperature and photope- 
riod influence the form of the phosphorus quota curve as well as 
the share of phosphorus uptaken by each species in competition 
experiments. 

Methods 

Ethics statement 

No specific permissions were required to isolate the diatom 
strains used in this study at the given location. All samples were 
taken from public waters. No field studies were performed 
involving endangered or protected species. 

Algal strains 

Growth experiments were performed under laboratory condi- 
tions for Stephanodiscus minutulus (Kiitz.) Cleve and MoUer 
(BaciUariophyceae), strain Mue0511A6, isolated in 2005 by 
Nicklisch, and Mtzschia acicularis W. Smith (BaciUariophyceae), 
strain Mue070319Cl, isolated by Nicklisch in 2007. All strains 
were isolated from Lake Mtiggelsee (Berlin, 52.44°N 13.65°E) and 
were axenic. 

Algae Cultivation 

S. minutulus and jV. acicularis were grown under P limitation in 
semi-continuous culture according to the chemostat principle [43] 
and under P-replete conditions according to the turbidostat 
principle. Algae were cultivated in FW04 medium, a fuUy synthetic 
freshwater nutrient solution according to Nicklisch et al. [1 1] with 
an ionic-composition similar to the water of Lake Mtiggelsee, 
supplemented with a trace element solution according to Nicklisch 
[43] and a vitamin solution according to GuiUard & Lorenzen [44] 
slightly modffied (final concentrations: 1 fjg L~' biotin, 1 ;Ug L ^ 
cobalamin and 100 //g L ' thiamine). The nutrients (300 /iM Si, 
200 /uM N, 10 l^M P, 2 Fe) were not limiting in the nutrient 
replete experiments at the low algal biomass densities we used. For 
P-limitation experiments, the P-concentration was decreased to 
1.2 fjM P. The solution was allowed to equilibrate with air by 
shaking to get a pH of about 8.3 at 20°C and then sterilised by 
filtering it through a membrane of 0.2 jUm pore diameter. The 
solution was subsequently heated by microwave to just below 
100°C but not allowed to boil. 

Algae were cultivated in Erlenmeyer flasks of 300 mL filled with 
100 mL suspension on an orbital shaker (50—90 rpm) in an 
incubator at 10, 15 and 20°C (±0.5°C). Irradiance was supplied 
from the side by fluorescent tubes of light colour Biolux and Warm 
White (Osram, Munich, Germany) at a ratio of 1:1 and cultures 
were positioned relative to the light source such that each culture 
received the same irradiance. The scalar photon flux density of the 
photosynthetically available radiation (PAR in fjmol quanta 
m ^ s ') was measured using a spherical sensor (QSL-101, 
Biospherical Instruments, California, USA). The light exposure 
in mol quanta m d was calculated as the sum of the irradiance 
over the photoperiod. Self-shading was minimised by the low 
biomass concentration (<300 /xg Chi a L ') and the shallow depth 
of the culture suspension in the flasks. 
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P-limited growth and competition experiments 

In P-limited growth experiments, cultures were grown at 10, 15 
and 20 °C with a photoperiod of 12 h d ' to investigate 
temperature interactions, and under photoperiods of 6, 9 and 
12 h d~' at 15°C to investigate the photoperiod interactions. In 
the text, we use the term irradiance to refer to the instantaneous 
photon flux density, and the term light exposure to refer to the 
total daily light supply, often called the daily dose. We avoid the 
use of dose because this is sensu stricto absorbed light [1 1,45]. Light 
exposures were saturating or near saturating for all experiments 
(all growth rates were S80% of the asymptotic maximum growth 
rates) as determined from growth-irradiance curves for these 
species measured at the experimental temperatures and photope- 
riods [11,36,38]. Irradiances in the experiments were adjusted 
to maintain a roughly constant daily light exposure of 3-5.6 mol 
m d . First, nutrient replete culture's were acclimated to the 
desired conditions of temperature, photoperiod and irradiance for 
5 to 7 days. During this phase and the following experiment, the 
biomass and growth rates were determined every day or every 
second day (at low growth rates) as described below. After 
acclimation when the cultures had reached balanced growth, the 
maximum specific growth rates were determined in turbidostat 
mode. Here cultures are diluted to a constant starting biomass 
each day and the dilution rate is adjusted accordingly. The growth 
rate is then determined from the dilution rate. The growth 
experiments were carried out for two to four weeks to get a mean 
of the specific growth rate. The longer measuring periods were 
needed for low growth rates at low temperatures or short 
photoperiods. 

To determine P-limited growth rates, nutrient replete cultures, 
which had been acclimated to the experimental conditions as 
described above, were diluted with P-free medium in such a 
proportion to achieve a final concentration of 1.2 //M P and 
subsequently allowed to starve for 2 to 5 days. The algae were then 
further cultivated semi-continuously with a nutrient solution 
containing 1.2 /iM P at a constant dilution rate (0.2-0.8 d~'). 
The cultures reached a steady state within 1 to 2 weeks so that the 
specific growth rate was equal to the dilution rate and biomass was 
constant at each dilution. It was assumed that the P added at 
dilution is absorbed very rapidly and that biovolume then 
increases practically linearly. The P-quota associated with the 
measured growth rate can then be calculated as the total P 
concentration per mean biovolume between dilutions. Here the 
actual short term growth rates measured between dilutions, which 
usually differed slightly from the applied dilution rates, were 
considered. The maximum yield (reciprocal of the minimum quota 
Qo) was determined by allowing the cultures to grow without 
dilution until they reached a maximum biomass (usually 1-2 
weeks). The total P-content of cultures was checked at the end of 
the experiments by standard methods [46]. 

In competition experiments, axenic, unialgal cultures of S. 
minutulus and jV. acicularis were first acclimated to the desired 
conditions and then grown under P-limitation at the dilution rate 
required for the later competition experiments. Once the cultures 
had reached balanced growth within 1-2 weeks, equal culture 
volumes of the two species were mixed so that each species initially 
contained 50% of the total phosphorus. Occasionally different 
proportions were used. The mixed cultures were further cultivated 
semi-continuously for one to three weeks at a constant dilution rate 
of 0.4-0.5 d"' at 10, 15 and 20°C with a 12 h d"' photoperiod 
and also at 15°C with a 6 h d~' photoperiod. The proportions of 
each species were measured in fixed samples (0.25% glutaralde- 
hyde final concentration) every one or two days using a flow 
cytometer (FACStar-Plus, Becton Dickinson, San Jose, USA) 



according to Nicklisch and Steinberg [47]. While there was some 
error in the absolute cell counts by the flow cytometer, the 
proportions of the species were accurate as verified by microscopic 
counts. Cells from the different species in the mixed cultures were 
distinguished based on their forward scatter signals and the Ghl a- 
autofluorescence excited by a blue laser at 488 nm. In competition 
experiments, the competitive ability (CA) of the weaker species {S. 
minutulus) relative to the stronger species (jV. aciculmis) was 
quantified by 

CA=1 + ^ (1) 

where /inet is the net specific growth rate (£0) and D is the applied 
dilution rate. The competitive ability is thus 1 when the two 
species coexist and 0 when the weaker species does not grow and is 
washed out at the dilution rate. 

Biomass and specific growth rate measurements 

During experiments, cultures were sampled every one or two 
days, always 3 hours from the beginning of the photoperiod and 
the biomass was determined by photometry at a wavelength of 
436 nm (5 cm cuvette, Shimadzu photometer type. UV-2401 PC). 
The focus of the light beam and the distance between cuvette and 
photomultiplier in this photometer excludes most of the scattered 
light, including forward scattered light, from detection. Therefore, 
the measured absorbance was caused by scattering (about 80%) 
and pigment absorption (about 20%). In parallel, chlorophyll 
fluorescence (F„) and variable fluorescence (F^) were measured 
using a Xenon-PAM Fluorometer (Heinz Walz GmbH, Effeltrich, 
Germany) after dark adaptation for at least 20 minutes. F^ is 
closely correlated to chlorophyll a content and F„ which is the 
increase in fluorescence above F^ after a light saturation pulse 
(Kautsky effect), is closely related to total PS II activity [48]. The 
chlorophyll and variable fluorescence were used to monitor the 
condition of the cultures and also as surrogates for biomass. 

Specific grovrth rates were thus calculated separately using the 
change in absorbance at 436 nm, and F^ The growth rates 
calculated with these three parameters are equal when the culture 
is in quasi-steady state. The three growth rate values (average 
difference 1 %) were averaged to obtain the most accurate estimate 
of the true rate. Absorbance was used rather than direct biovolume 
measurements due to higher repeatability and shorter measuring 
and dilution times, which minimised culture disturbance. Before 
taking samples, we performed absorbance and fluorescence 
measurements to ensure that growth was balanced, in which case 
cultures are in quasi steady state. Here the biomass composition is 
nearly constant and the specific growth rates are not only biomass- 
specrfic but also C- or Chi a-specific, provided that the cultures are 
measured at the same time within the photoperiod. 

Phosphorus, carbon and chlorophyll content 

The phosphorus content of cells was determined as the total 
phosphorus concentration in the fresh medium divided by the 
biovolume concentration of the culture at sampling, assuming that, 
under P-limitation and axenic conditions, all phosphorus was 
incorporated into the cells. Biovolume was measured with a cell 
counter (CASY, Modell TTC, Scharfe System, Germany). Carbon 
c()nt(;nt was measured with a C/N analyser (Vario EL, Elementar 
Analysensysteme, Hanau, Germany) after collecting cells on pre- 
rinsed, pre-fired and pre-weighed filters (Munktell 25 mm). 
Chlorophyll a was measured for each culture at each dilution 
rate and experimental treatment with a Waters Alliance (MUford, 
USA) high performance liquid chromatography (HPLC) system 
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with photodiode array detector, according to the method of 
Shatwell at al. [38]. In addition, each time a culture was sampled, 
chlorophyll content was measured using the chlorophyll fluores- 
cence, F„, calibrated against the HPLC measurements. Chloro- 
phyll content for each culture was then taken as the mean of the 
HPLC and fluorescence measurements. 

Statistics and model fitting 

To analyse P-limited growth measurements, we fitted and 
compared different types of rectangular hyperbolic (HYP) and 
exponential (EXP) saturation curves, which describe the specific 
growth rate (ji) under P-limitation as a function of P-quota (Q): 



Qm(Q-Qo) 
'Q(Qm-Qo) 



HYP, after Droop [26] 



(2) 



l-exp(-ln2^|-l 
l-exp( -ln2(^-l 



EXP, after Fuhs [49] (3) 



Here Qj) is the subsistence or minimum cell quota. The original 
Droop and Fuhs models were formulated in terms of the 
theoretical maximum growth rate at infinite quota (jl'm) but here 
have been reformulated using the real maximum specific growth 
rate {/v) which occurs at the maximum quota (Qm)- 

In addition, we also fitted normalised quota models which 
contain a dimensionless shape parameter (A"Q, k) in place of a half 
saturation constant: 



(i+^e)(Q~Qo) 



'Q-Qo+^e(Qm-Qo) (4) 
HYP, after Flynn [50] based on Caperon & Meyer [51] 



l-exp(-.^|-l 



1 — exp 



EXP, based on Fuh [49] (5) 



The normalised shape parameters are related to the half- 
saturation constant (kg) as follows: 



Eq. 3. KQ^and k are not direcdy comparable with each other and 
they should not be confused with Icq, which has the same units as 
Q. KQ^and k do however provide a basis for generalising the form 
of the rjuota curve irrespective of the units of biomass used [27]. 
Decreasing KQ^ and increasing K (at constant fl^) indicate an 
increase in efficiency of P usage. 

To analyse the effect of temperature and photoperiod on P- 
limited growth, we fitted Eq. 5 to each experimental treatment and 
examined the effect of temperature and photoperiod on the fitted 
model parameters. To determine whether the temperature or 
photoperiod effect was significant, we first fitted Eq. 5 to the 
pooled data from all temperature (or photoperiod) treatments. We 
then added a linear temperature (or photoperiod) dependence to 
the model parameters and refitted the model. These more 
complex models were subsequently compared with the simpler 
model using F-tests according to Bates and Watts [52]. If the 
residual variance of the more complex model with the temperature 
(or photoperiod) dependence was significantly lower than the 
residual variance of the simpler model without a temperature (or 
photoperiod) dependence (after accounting for the increased 
model complexity), we inferred that the temperature (or photo- 
period) effects were significant. Eqs 2-5 were compared with 
Akaike's Information Criterion (AIC), which was interpreted 
according to Bumham and Anderson [53]. Although we fitted 
models using both cell P-quota and biovolume P-quota, we only 
show results for biovolume quota because the fits were similar and 
the biovolume quota allows better comparisons between the 
sp(-{ i(;s. Statistical analyses, curve fitting and model simulations 
(described below) were performed using R version 3.0.1 [54] with 
the packages deSolve [55] and FME [56]. 

Model and simulations of relative P uptake rates 

To estimate the P- uptake rates of 5. minutulus and jV. acicularis, we 
simulated the growth of these species in the competition 
experiments described above in semi-continuous P-limited culture. 
The model parameters are described in Table 1. The simulation 
was run in steps, from one dilution to the next. Between dilutions 
the two species grew as in a batch culture, such that 



dX 
dT 



dQ 
dt 



= HX 



(8) 



(9) 



KQ- 



Qm-Qo 



(c= ln2 



Qo 



(6) 



(7) 



Since the Droop and Fuhs models are half saturated at 2Qo (i.e. 
kc) = Q)), the Droop relation applies when KQj= Qs^/iQjn - Qji) 
(Eq. 6) or when K = In 2 (Eq. 7). Accordingly in these cases, Eq. 4 
becomes equivalent to Eq. 2 and Eq. 5 becomes equivalent to 



Because the maximal uptake rate (V^) is much higher than 
required to satisfy immediate needs for growth if cells are 
substantially P-limited [15,17], we assumed the P added during 
dilution was absorbed instantly and completely, so that the 
dissolved P concentration in the medium was negligible. Further- 
more, we assumed that the added P was distributed between the 
two species (i and j) in fixed proportions, Vrei, i and Vrei, j, such that 
Vrei, i+Vrei_ j = 1 . Therefore, at the end of each step the biomass (X) 
and P-quota (Q) of each species were modified to account for 
dilution as follows: 



X*=/X 



(10) 



PLOS ONE I www.plosone.org 



4 



July 2014 I Volume 9 | Issue 7 | e102367 



Factor Interactions in Diatoms under P-Limitation 



/ = e-°^' (12) 

where X* and Qf denote the new biomass and P-quota, 
respectively, at the start of the next simulation step directly 
following dilution, S^^ is the P concentration in the fresh medium, 
and /is the dilution factor, which depends on the dilution rate D 
and the time until the next dilution. At (Eq. 13). Note the 
subscripts, i and j, which have been introduced because the 
nutrient absorbed by one species depends on the biomass of the 
other. 

The variables V„:i_ i and V„.i j are biomass-specific and are 
therefore proportional to the P-uptake rates of the two species 
relative to each other. In effect, our approach assumes a linear 
relationship between external nutrient concentration and uptake 
rate. The slope of this linear relationship (Vrd, i and V^ei^ j) is 
proportional to the uptake affmity, or the initial slope of the 
Michaelis-Menten uptake curve V^/k^, where k^, is the half- 
saturation coefficient [42,57]. In a calibration routine during the 
simulations, i and V^ei, j were fitted to the measured data in 
each competition experiment to investigate how the relative 
uptake rates change with temperature and photoperiod. 

Results 

P-replete growth 

The different models (Eqs. 2-5) were fitted to the growth rates 
of jV. ackularis and S. minutulus under P-limited and P-replete 
conditions. Eq. 5 performed the best overall for both species and 
the Droop model (Eq. 2) performed worst (Table 2). The 
normalised, three-parameter models (Eqs. 4, 5) performed better 
for JV. ackularis whereas the exponential models (Eqs. 3, 5) tended 
to perform better for S. minutulus. Therefore, the normalised 
exponential model Eq. 5 was used to compare temperature and 
photoperiod treatments for the following analysis (Figure 1). 



The maximum (P-replete) specific growth rates increased with 
increasing temperature and photoperiod for both species as 
expected (Figure 2). M. ackularis had a higher maximum specific 
growth rate (JJ^) than S. minutulus at 10-20 °C and 12 h d ' 
photoperiod but growth rates were similar under shorter 
photoperiods. 

The maximum P-quota (Q^J when cells were growing at their 
maximum rate (//„,) under luxury consumption decreased with 
increasing temperature from about 5 //g P mm ^ at 10°C to 3 /Jg 
P mm ^ at 15 °C and 20 °C for S. minutulus (p<0.001. Table 3). 
Overall JV. acicularis had a higher P-quota than S. minutulus (p< 
0.001), with Qp, decreasing significantly (p<0.001) from 6 ^Ug P 
mm"'^ at 10 °C to 4 ^g P mm"^ at 20 °C. 

Temperature and photoperiod effects on P-limited 
growth 

The minimum P-quota required for growth (Qp) was lower for 
S. minutulus than for jV. ackularis, and the normalised shape 
parameter K, which provides an indication of the efficiency of P 
usage, was higher for S. minutulus (Figure 2C, H). This shows that S. 
minutulus is more efficient and can produce more biomass from a 
given amount of internal P than jV. ackularis. Qj) was independent 
of both temperature and photoperiod for both species (p&0.39, all 
cases). Although model fits suggested that was significantiy 
higher for .S". minutulus at 20°C than in the other experiments 
(Table 3), there was no significant difference between temperature 
treatments when the quotas measured in cultures grown to 
stationarity (u = 0) were compared direcdy by ANOVA (p = 0.43). 

The normalised shape parameter K did not change significandy 
with photoperiod for either S. minutulus (p = 0.80) or jV. acicularis 
(p = 0.09, Figure 2H). k increased significantly with increasing 
temperature in S. minutulus and jV. acicularis (p^O.OOl in both cases, 
Figure 2C). We also checked the temperature and photoperiod 
dependence of the normalised parameter ^Q,in Eq. 4, which has a 
somewhat different meaning to K because it includes the 
temperature dependent Q^^. KQ_ decreased with increasing 
temperature (p = 0.03) and increased with photoperiod (p = 0.02) 
for jV. acicularis but was constant over temperature and photope- 
riod for S. minutulus (Figure 2D, I). The half saturation coefficient 
kq, which can be calculated from Eq. 7, overall decreased 
significandy with increasing temperature in both species (p£ 
0.001), although in S. minutulus, kq^was not significantly different at 
15 and 20°C (p = 0.19. Figure 2E). The temperature dependence 
of kqwas substantially weaker (flatter in Figure 2E) for S. minutulus 



Table 1. Model parameters and variables for simulation of P-competition and relative uptake rates in semi-continuous culture 
experiments. 





Parameter 


Meaning 


Units 


X 


Biovolume 


mm^ L"' 


Q 


P-quota 


ftg P mm"^ 




Concentration of P in fresh medium 




Vrel, 1 


Biovolume-specific proportion of P absorbed by species i 


Dimensionless 




Specific growth rate, given by Eq. 5 with the parameters in Table 3. 


d-^ 


F 


Dilution factor (proportion of culture retained at dilution) 


Dimensionless 


T 


Time 


d 


At 


Time until next dilution 


d 


D 


Dilution rate 


d-' 


Subscripts / and j refer to the two species. 
doi:l 0.1 371/journal.pone.Ol 02367.t001 
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Stephanodiscus minutulus Nitzschia acicularis 



T ' 1 ' 1 ' 1 ' 1 ' 1 ' p ' 1 ' 1 ' 1 ' r 




Figure 1. Specific growth rates as a function of biovolume P-quota at different temperatures and photoperiods. Points show 
measured values for S. minutulus (A-E) and N. acicularis (F-J), solid lines show the fitted model (Eq. 5) with the parameters in Table 3. The dotted lines 
show the maximum specific growth rates (fi^) predicted by the model of Nicklisch ef al. [11] at the respective temperature, photoperiod and light 
exposure. 

doi:1 0.1 371/journal.pone.01 02367.g001 



Table 2. Comparison of different models fitted to experimental results. 





Experimental conditions 


AAIC 








Temperature (°C) 


Photoperiod (h d ^ 


Eq. 2 after Droop [26] 


Eq. 3 after Fuhs [49] 


Eq. 4 after Flynn [50] 


Eq.5 


Stephanodiscus minutulus 


10 


12 


41.4 


39.1 


0 


2.7 


15 


12 


71.4 


16.0 


7.3 


0 


20 


12 


10.7 


0 


5.2 


2.0 


15 


6 


46.2 


2.2 


15.8 


0 


15 


9 


20.7 


0 


33 


1.1 


Nitzschia acicularis 


10 


12 


80.9 


75.3 


0.2 


0 


15 


12 


91.3 


85.3 


3.5 


0 


20 


12 


19.2 


13.3 


0 


2.7 


15 


6 


36.9 


17.4 


5.5 


0 


15 


9 


44.8 


37.0 


0 


3.9 



Lower AAIC values indicate better model fits {best model = 0). AAIC <2 indicates substantial support that the model in question is the best model; AAIC >1 0 indicates 

essentially no support [53]. 

doirl 0.1 371 /journal.pone.Ol 02367.t002 
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Figure 2. Temperature and photoperiod dependence of physiological growth parameters for 5. minutulus and N. acicularis. The 

temperature dependence (A-E) was measured at 12 h d"^ photoperiod and the photoperiod dependence (F-J) at 15 C. The physiological 
parameters maximum specific growth rate (1;^), minimum P-quota (Qo), normalised shape parameter (k) and half saturation coefficient (kg) were 
fitted using Eq. 5, and the normalised shape parameter [KQ) was fitted using Eq. 4. 
doi:1 0.1 371 /journal.pone.01 02367.g002 



than jV. acicularis. In S. minutulus, K was approximately equal to In 2 
{ = 0.69) (except at 1 0°C, where k = 0.32), whereas in JV. acicularis k 
was lower than In 2 by a factor of 2 to 4 (Table 3). Therefore the 
simpler Droop (Eq. 2) and Fuhs (Eq. 3) models, which are half 
saturated at 2Qj) (i.e. such that K = ln 2), provided a reasonably 
good description of the growth of S. minutulus but not of jV. 
acicularis. 

Model of factor interactions 

To describe the interactions between temperature, photoperiod 
and P-limited growth, we used the model of Nicklisch et al. [11], 
which considers the interactions between light exposure, photo- 
period and temperature, to calculate the maximum growth rates 
{]lra), and coupled it to a normalised quota model to account for P- 
limitation. The maximum growth rates predicted by the model of 
Nicklisch et al. agreed very well with the maximum growth rates 
measured under P-replete conditions (dotted line in Figure 1) for 
both species. For the quota model, the temperature and 
photoperiod dependencies of the model parameters shown in 



Figure 2 were formulated in model terms. Because Qj) was 
independent of photoperiod and temperature in both species, it 
was assumed constant in the quota model formulations. Therefore 
we tested Eq. 4 with constant and Eq. 5 with temperature 
dependent k as follows: 

'v = /v'opf/(T) (13) 



Where K^pt is the normalised shape parameter at optimum 
temperature and7(T) is the temperature function after Lehmann 
et al. [58] used in the model of Nicklisch et al. [1 1] with 
temperature T, optimum temperature T„pt and minimum 
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temperature T^i^ (aU °C). Fits of both normalised quota models to 
the whole data set for each species showed that there was very 
strong evidence (AAIC >10) that model 5 with a temperature 
dependent k {S. minutulus: AIC=— 822.1, M ackularis: AIC = — 
358.4) described the measured growth rates of each species better 
than model 4 with constant Kd{S. minutulus: AIC = -809.8, N. 
acicularis: AIC = —328.8). This is in agreement with the model 
comparisons of the individual experiments described above, where 
model 5 produced the best fit overall. 

Therefore we adopted model 5 with the parameters in Table 4, 
which overall provided a good fit, and used this model to plot the 
response surfaces under combined limitation by phosphorus and 
temperature or photoperiod (Figure 3). S. minutulus could achieve 
the same growth rate as jV. acicularis at considerably lower 
temperatures, photoperiod and P-quota combinations. A horizon- 
tal cross section through the response surface shows growth 
isoclines (constant growth rate) as in TUman's [59] interaction 
diagrams. The curved regions of these growth isoclines indicate 
that an increase in photoperiod and temperature can compensate 
a decrease in P-quota and the regions of interaction correspond to 
the prevailing spring conditions for these species (e.g. roughly 8- 
15°C and 4—8 h d~' photoperiods at a growth rate n = Q.b). 

Competition experiments and relative P uptake rates 

In competition experiments under P-limitation, N. acicularis was 
the stronger competitor at all temperatures and photoperiods 
tested. The experiments showed that, relative to jV. acicularis, the 
competitive ability of S. minutulus increased with decreasing 
temperatures under a 12 h d~' photoperiod (Table 5, p<0.001, 
multiple linear regression on time and temperature) and was 
higher at 6 h d~' photoperiod than at 12 h d~' photoperiod at 
15°C (p<0.001, t-test on linear regression slopes). 

The relative uptake rates were estimated from the model 
simulations, which accurately described the measured cell 
proportions of the two species in the competition experiments 
(Figure 4). The estimated relative uptake rates (Figure 5) showed a 
similar pattern to the competitive ability results shown in Table 5. 
The simulations showed that the uptake rate of S. minutulus relative 
to M. acicularis was higher at 10°C compared to 15 and 20°C 
(Table 5, p = 0.0004, ANOVA). At 10°C under a photoperiod of 
12 h d ', S. minutulus absorbed about 25-30% of the added P, 
whereas at 20°C it absorbed around 10%, while jV. acicularis 
absorbed 90%. There was no significant difference between the 
relative uptake rates at 6 h d~' and 12 h d~' photoperiod 
(p = 0.6), where S. minutulus and JV. aeicularis absorbed on average 
14% and 86% of added P, respectively. 

Carbon content 

The mean carbon contents of jV. acicularis and S. minutulus were 
148±39 /zg C mm-^* (n= 15) and 132±18 /ig C mm"^ (n= 13), 
respectively. In S. minutulus, the carbon content per unit biovolume 
did not change with dilution rate under P-limitation (p = 0.28). In 
jVl acicularis carbon content per unit biovolume increased with 
increasing dilution rate under P-limitation according to 
C:X= 1 18xD+56 (R^ = 0.82, p<0.01), where C:X is die carbon: 
biovolume ratio {jig C mm ^) and D is the dilution rate (d '). 

Clilorophyll a content 

Chlorophyll a content increased slightly with increasing dilution 
rate (decreasing P limitation) up to 0.7 d in S. minutulus and was 
relatively constant over this range for N. acicularis (Figure 6). 
Chlorophyll a was higher in P replete cultures (dilution rate above 
0.7 d ') than in P limited cultures in both species. The mean 
chlorophyll a content of S. minutulus and N. acicularis under P 



limitation was 2.38±0.77 ng mm {n= 195) and 4.33±1.13 ng 
mm~^ (k=115), respectively. Under nutrient replete conditions, 
the chlorophyll a content was 4.59±1.36 /ig mm (n = 28) and 
6.28 ± 2.49 //g mm ^ (h=31), respectively for the two species. 
Under P limitation, chlorophyll a did not significantly depend on 
temperature for either species (p>0.28) but did depend on 
photoperiod for both species (/><0.001), where the content was 
higher under short photoperiods than under long photoperiods. 

Discussion 

Our study demonstrated that temperature influenced the 
relationship between growth and P-quota of S. minutulus and M. 
acicularis, but the effect of photoperiod length at constant daily light 
exposure was not significant. Moreover the effect of temperature 
differed between species in a way that a fixed quota curve like the 
Droop model could not entirely account for. Similarly temperature 
also influenced the relative phosphorus uptake rates of the two 
species but the photoperiod apparently had no influence here. We 
therefore confirm our hypothesis that there is a complex 
interaction between temperature and phosphorus limitation but 
reject our claim that the photoperiod has an effect. 

P-replete growth and model comparison 

In our growth experiments, the measured P-replete growth rates 
(jj^ increased non-Unearly as expected with increasing tempera- 
ture and photoperiod. These growth rates agreed very well with 
the predictions of a previously published model [11]. 

One reason for the Droop model's poor performance compared 
to the other tested models is that the fixed form of the Droop 
equation could not adequately describe the growth rates of N. 
acicularis because K:<ln 2. There are other examples where the 
Droop equation did not apply as evident from non-linear 
relationships between n and 1/Q, [9,30,32,60,61], which typically 
occur when the growth-quota cur\ (- has a flatter form such that 
K<ln 2. This explains why the 3-parameter models (Eqs. 4, 5) 
fitted better than the 2-parameter models (Eqs. 2, 3) for jV. 
acicularis. Although the 3-parameter (variable form) curves also 
provided a significantly better fit for S. minutulus than the 2- 
parameter (fixed form) curves, a 2-parameter model would in our 
opinion be adequate for this species. 

We used biovolume to represent biomass as the basis for P- 
quota estimates. In diatoms, which can be quite vacuolated, 
biovolume can lead to biased estimates of the biomass. An 
alternative is to use C as the quota basis, which is useful as a 
common biomass currency in modeUing studies. However, C- 
content can also lead to biased estimates of biomass under nutrient 
limitation when excess photosynthates are stored, such as 
chrysolaminaran in diatoms [62] . Dry weight poses difficulties in 
diatoms due to their silicate frustules and ceU number can cause 
problems due to variability of cell size [27]. We chose biovolume 
because it is fast and easy to measure and it allows modelling 
estimates of individual species to be verified with field samples, 
which is not possible with C based quotas. In our cxpcrimc-nts, the 
results were similar using biovolume based and cell based P- 
quotas, and for S. minutulus also C based quotas because C content 
per unit biovolume was constant. In ,¥. acicularis, C per biovolume 
increased with dilution rate. Using a C based P-quota would 
therefore result in a more linear quota curve and further decrease 
K below In 2. 

Temperature 

Our results showed that the P-quota increased with decreasing 
temperature for a given growth rate as evident from a 
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A S. minutulus C S. minutulus 




B N. acicularis D N. acicularis 




Figure 3. Interaction diagrams sKiowing interactions between P-quota, temperature and photoperiod. Diagrams sinow the response of 
specific growth rate (d^^) to the combined effect of temperature (°C) and P-quota (fjg P mm"^) at 12 h d^' photoperiod for S. minutulus (A) and N. 
acicularis (B), and also to the combined effect of photoperiod (h d"') and P-quota at 15'C for S. minutulus (C) and N. acicularis (D). 
doi:1 0.1 371 /journal.pone.01 02367.g003 



Table 4. Model of factor interactions (Eqs. 5, 14, 15) and corresponding parameters for P-limlted growth of S. minutulus and N. 
acicularis. 





Parameter 


Description 


Units 


5, minutulus 


N. acicularis 




Maximum (P-replete) specific growth rate 


d-^ 


a 


a 


Qo 


Minimum P-quota 


/ig P mm"^ 


0.43 (0.41-0.44) 


0.53 (0.49-0.57) 


Kopt 


Normalised shape parameter at optimum temperature 


Dimensionless 


0.69 (0.64-0.74) 


0.34 (0.29-0.39) 


Topt 


Optimum temperature 


T 


20.7° 


21.7' 


"l"mln 


Minimum temperature 


°C 


-o.e 


1.0' 


RSE 




d-' 


0.070 


0.092 


df 






330 


185 


^Calculated or adopted from Nicklisch ef ol. [1 1]. 

RSE: residual standard error, df: degrees of freedom, 95% confidence intervals shown 
doi:1 0.1 371/journal.pone.01 02367.t004 


in parentheses. 
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Figure 4. Simulation of P-competition experiments between N. 

acicularis and S. minutulus. This example experiment was performed 
at 15 C, 12 h d^^ photoperiod and 3.7 mol PAR m"'^ d^V Lines show 
the simulated values of (A) proportions of each species by biovolume, 
(B) the biovolume P-quotas, and (C) specific growth rates. Points show 
measured values. N. acicularis was the stronger competitor (solid line) 
and S. minutulus the weaker competitor (dashed line). 
doi:10.1 371/journal.pone.01 02367.g004 

temperature-dependent k and Q,,,, which is in accordance with 
other studies [13,33]. The normahsed sliape parameter k clearly 
increased with decreasing temperature in both species, which 
indicates that the shape of the quota curve changes with 
temperature. We showed that the temperature dependence of k 
is approximately the same as the temperature dependence of 
nutrient replete growth (c.f Eqs. 14, 15), although it is unlikely to 
be exactly the same [30] . In absolute terms the increase in with 
decreasing temperature was only small (but nevertheless signifi- 
cant) for S. minutulus. An increase in k^ corresponds to a decrease 
in growth rate at a certain P-quota, which suggests that the much 
smaller temperature effect on k^ may reflect an adaptive strategy 
of S. minututus because it is a cold-adapted, early spring species 
whereas jV. acicularis is adapted to warmer temperatures and longer 
photoperiods and typically grows in late spring [11,38]. 

The absence of a temperature dependency of the minimum 
quota Qj) in our study is not consistent with other findings that Qp 
increases with decreasing temperature [13,30,32], although there 
are some examples where Qj) did not increase with decreasing 
temperature [33,34,49]. A decrease in cell size with increasing 
temperature [13], which we also observed in our study (Table 3), 
might partially explain this discrepancy if nutrient quotas are given 
on a cell basis [14], but there does not seem to be a consistent 
pattern in the literature to provide a clear answer to this. Another 
explanation could be that the discrepancy is due to methodological 
differences or curve fitting because the choice of model 
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Figure 5. P-uptake rates of 5. minutulus relative to N. acicularis. 

The relative uptake rates were measured at different temperatures with 
a 12 h d"' photoperiod (A) and at different photoperiods at 15 C (B). 
Open circles show the fitted values of parameter Vrei for S. minutulus, 
error bars show 95% confidence intervals of the fitted values. The 
dimensionless relative uptake rates (Vrei) shown refer to the proportion 
(%) of added P absorbed by S. minutulus. It follows that the proportion 
absorbed by N. acicularis is 1- V,ei of S. minutulus. 
doi:1 0.1 371/journal.pone.01 02367.g005 



considerably affects the parameter values. Furthermore, if the 
Droop equation is fitted, then a change in Qj) with temperature 
win also reflect a change in K or Icq due to the fixed form of the 
quota curve, so that a temperature dependence of Qp found by 
model fitting could be just a model artefact. It is therefore 
interesting to note that studies that found a temperature 
dependence of Qj) estimated this as the axis intercept extrapolated 
from a fitted curve [13,30,32]. On the other hand, the studies that 
found no temperature dependence measured Qj) directly from the 
maximum yield of cultures grown to stationary phase [33, 34; this 
study] . One exception was the study of Fuhs [49] , who estimated 
Qp by extrapolation but nevertheless did not find a temperature 
dependence. Ahlgren [14] also questioned the significance of the 
temperature dependence of Qp for phosphorus in some studies. 
We checked the possibility that the temperature dependence of Qp 
could be an artefact of model fitting by removing the data points 
measured at statioiiarity ljil — 0) and refitting the Droop equation to 
our data. This resulted in a highly significant temperature 
dependence of Qp for jV. acicularis but no temperature dependence 
for S. minutulus. This is consistent with the need for additional P- 
rich ribosomes [63] to compensate a decrease in temperature and 
still maintain the same protein synthesis and growth rate 
[13,47,64]. However, this may not apply when the growth rate 
is zero. We suggest that further research is necessary to clarify 
whether Qp is temperature dependent because this is relevant for 
optimum N:P ratios and overall phytoplankton stoichiometry [65] . 

The majority of published results (e.g. [13,32]) and our results 
for JV. acicularis indicate that the form of the quota curve changes 
with temperature, making a fixed-form curve such as the Droop 
model inappropriate in most cases to account for temperature 
interactions with phosphorus [30]. The fact that the Droop 
relation applied reasonably well for S. minutulus and the 
temperature interactions with phosphorus-limited growth could 
be relatively well described by a multiplicative model with a fixed 
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Figure 6. ChloropKiyll a content at different levels of P limitation (dilution rate). Different symbols show temperature treatments at 
12 h d"' photoperiod (A, B) and photoperiod treatments at 15°C (C, D) for S. minutulus (left panels) and N. acicularis (right panels). Symbols have 
been staggered horizontally to avoid overlap; dilution rates increased from 0.2 d"^ to 0.7 d"' at 0.1 d"^ intervals. Higher dilution rates show 
maximum growth rates. Error bars ±1 SD. 
doi:1 0.1 371/journal.pone.01 02367.g006 



PLOS ONE I www.plosone.org 



12 



July 2014 I Volume 9 | Issue 7 | e102367 



Factor Interactions in Diatoms under P-Limitation 



curve is presumably due to specific adaptation of this cold-water 
species and serves to higlilight the species-specific nature of the 
effect of temperature on P-limited growth. 

Photoperiod 

The lack of a significant effect of the photoperiod on Qji or K 
indicates that the photoperiod does not alter the shape of the quota 
curve for S. minutulus or M. acicularis. Further, it suggests that the 
effects of photoperiod and P-limitation are multiplicative and that 
a decrease in the photoperiod would require an increase in P- 
quota to maintain the same growth rate, as was evident from the 
interaction plots (Figure 3). This is consistent with the results of 
Litchman et al. [20], who found that the response of several 
phytoplankton species to daylength depended on the P-status and 
vice versa. A similar result was also found for Limnothrix (formerly 
Oscillatoria) redekei because and Qp were the same under both 
continuous light and a photoperiod of 1 2 h d~ ' [33] . On the other 
hand, Riegman and Mur [21] found only a narrow range of 
interaction between internal phosphorus and the photoperiod for 
Planktothrix agardhii. These authors suggested that Liebig's Mini- 
mum Law applies to P. agardhii because it is a shade adapted 
species and they hypothesised that high light adapted species 
might have a broader interaction range. We used our model, 
which contains a multiplicative interaction between P-quota and 
photoperiod, to examine the interactions of S. minutulus and jV. 
acicularis (which are slightly less shade-adapted than P. agardhii) 
under the same conditions as used by Riegman and Mur (results 
not shown). The interaction range was somewhat broader, 
supporting Riegman and Mur's hypothesis, but also demonstrating 
that a multiplicative interaction can produce visually similar results 
at the low growth rates these authors used (0.01 h '). Evidence 
shows that a decrease in irradiance (//mol photons m ^ s ') 
generally requires an increase in P-quota to maintain the same 
growth rate, suggesting an interaction between light and P- 
limitation [9,14]. Furthermore, Healey's [9] results demonstrated 
that Qp and for phosphorus (both cell- and C-based quotas) 
were independent of irradiance, indicating that the interaction 
between irradiance and phosphorus is multiplicative. This might 
help to explain the similar interaction with the photoperiod in our 
study, where k^ and K were independent of the photoperiod and, 
according to Healey's results, also the irradiance. Further evidence 
is provided by the fact that in S. minutulus and M. acicularis is 
unaffected by either irradiance or photoperiod (Figure 7) [66] . 

We found that chlorophyll content decreased under P 
limitation, which is consistent with other studies [31,67,68]. Since 
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Figure 7. Maximum P-quotas (Q,^) at different pKiotoperiods. 

The data shown are for Stephanodiscus minutulus and Nitzschia 
acicularis measured at 20 C in semi-continuous culture in P-replete 
medium, and were published in Giersdorf [66]. 
doi:1 0.1 371/journal.pone.01 02367.g007 



chlorophyll increases under light limiting conditions, photoaccli- 
mation as well as light-phosphorus interactions can play a role 
under combined nutrient and light limitation. This should be more 
relevant in eutrophic lakes, which through self-shading tend to be 
more turbid with a shallower euphotic depth and shorter effective 
photoperiod than in clear, strongly P limited oligotrophic lakes. 
Therefore, the type of interaction between phosphorus and 
photoperiod or light should be accounted for when coupling a 
quota model to a model of hght limitation. We do not expect that 
photoacclimation influenced our results because the cultures were 
grown at steady state. 

Competition and P-uptake 

To estimate relative P uptake rates, we employed a new 
approach which uses the relationship between growth and P-quota 
and the outcome of competition experiments. This approach 
allowed us to quantify the contribution of P-uptake and efficiency 
of P-usage to competitive abihty, but it could not produce absolute 
values of uptake rates. Our assumption of a simple linear 
relationship between uptake rate and external concentration can 
provide reliable results [69,70] and applies if the external P 
concentration is lower than the half-saturation coefficient of 
uptake (kjn). The P-concentration in our fresh medium (1.2 /iM P) 
was lower than measured k„-i values for most species provided by 
Gotham and Rhee [71] (mean 1.5 flM P, n = 6) and the external 
concentrations in our cultures can therefore generally be expected 
to be much lower than k„-,. Our linear relationship also ignores the 
feedback of the quota, which causes Vm to decrease as Q increases. 
However, the P-quota was always lower than Q^^ in model 
simulations for both species, except very occasionally after cultures 
were strongly diluted. Furthermore, this feedback depends on 
nutritional history rather than the instantaneous P-status [69,72] 
allowing phytoplankton to overshoot their maximal quotas in the 
short term [17]. Therefore, we assume that ignoring the feedback 
from the quota did not affect our results to any great degree. 

The outcome of the experiments showed that jV. acicularis was a 
strong competitor under P-limitation, which is consistent with 
other studies [43,73]. Interestingly S. minutulus, which was the 
weaker competitor under all experimental conditions, could use 
internal P more efficiently than jV. acicularis under most conditions, 
due to its lower Qq and higher K. The competitive advantage of jV. 
acicularis was due to a higher uptake affmity. Other studies found 
similar results, where the stronger of two species competing under 
nutrient limitation had higher uptake rates and a higher 
whereas the weaker competitor produced a higher yield from a 
given amount of nutrient [69,74]. A high maximum quota in 
combination with a high uptake rate enables a species to maintain 
a relatively high quota so that efficient use of internal P for growth 
is less important. This may suggest a trade-off between resource 
gathering and resource usage in biomass assembly [28] and 
stresses the importance of linking uptake and the quota curve [27]. 
It thus seems that the parameters of the quota curve alone provide 
litde information on the competitive abihty of a species without 
knowledge of the uptake kinetics. For this reason, Flynn [50] 
concluded that, in addition to the growth-quota relationship, 
models need to account for surge nutrient transport to properly 
describe competition between species. 

The influence of temperature on relative P-uptake rates was 
clear in the competition experiments. It was not possible to 
estimate the absolute temperature dependencies of uptake for each 
species, but it was possible to deduce how P was distributed 
between the species and how temperature and photoperiod 
influenced this distribution. S. minutulus increased its share of 
absorbed P from around 10% at 20°C to 25-30% at 10°C, 
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indicating that the uptake kinetics of S. minutulus are more cold- 
adapted than those of jV. aciculmis, which parallels the interaction 
between temperature and the growth-quota curve described 

above. 

The photoperiod appeared to have no influence on the relative 
uptake rates of the two test species. We cannot rule out that the 
photoperiod did alfect the absolute uptake rates of each species but 
that these effects cancelled each other out; however, we assume 
this is unlikely because there was no evidence that the photoperiod 
influenced other P-limited growth parameters like K, kfj or Qp. 
Sommer [75] found that the total daily light exposure rather than 
photoperiod or irradiance alone influenced the outcome of 
competition between marine phytoplankton under Si or N 
limitation. In our competition experiments, we tried to keep the 
daily light exposure constant but stiU found that the photoperiod 
influenced the outcome. However, this was obviously due to the 
effect of the photoperiod on maximum growth rates because 
relative P-uptake rates were unaffected. 

Altogether our study showed that temperature interactions with 
P-limited growth are complex and reflect species-specific niche 

References 

1. Jeppesen E, Sondergaard M, Jensen JP, Havens KE, Anneville O, et al. (2005) 
Lake responses to redueed nutrient loading - an analysis of contemporary long- 
term data from 35 ease studies. Freshwater Biol 50: 1747—1771. 

2. ShatwcU T, Kohler J, Nieklisch A (2008) Warming promotes cold-adapted 
phytoplanklou in temperate lakes and opens a loophole for OsciUatoriales in 
sprmg. (;iobal Change Biology 14: 2194-2200. 

3. Kirillin G (2010) Modeling the impact of global warming on water temperature 
and seasonal mixing regimes in small temperate lakes. Boreal Environ Res 15: 
279-293. 

4. Rinke K, Yeates P, Rothhaupt KO (2010) A simulation study of the feedback of 
phytoplankton on thermal structure via light extinction. Freshwater Biol 55: 
1674-1693. 

5. Jeppesen E, Moss B, Bennion H, Carvalho L, DeMeestcr L, et al. (2010) 
Interaction of climate change and eutrophication. In: Kernan M, Battarbee R, 
Moss B, editors. Climate Change Impacts on Freshwater Ecosystems. Oxford, 
UK: VVilcy-Blaekwell. 119-151.^ 

6. Kohlcr J, Hilt S, Adrian R, Nieklisch A, Kozerski HP, et al. (2005) Long-term 
response of a shallow, moderately Hushed lake to reduced external phosphorus 
and nitrogen loading. Freshwater Biol 50: 1639—1650. 

7. KiUiEim P, Kilham SS, Hecky RE (1986) Hypothesized resource relationships 
among African planktonic diatoms. Limnol Oceanogr 31: 1169-1181. 

8. Thompson P (1999) Response of growth and biochemical composition to 
variations in daylength, temperature, and irradiance in the marine diatom 
Thalassiosira pseudonana (Bacillariophyceae). J Phycol 35: 1215—1223. 

9. Hcaley FP (1985) Interacting effects of light and nutrient limitation on the 
growth rate oi Synechococcus linearis (Cvanophyeeae). J Phycol 21: 134—146. 

10. Foy RH, Gibson CE (1993) The influence of irradiance, photoperiod and 
temperature on the growth kinetics of three planktonic diatoms. EurJ Phycol 28: 
203-212. 

11. Nieklisch A, ShatweU T, Kohler J (2008) Analysis and modelling of the 
interactive effects of temperature and light on phytoplankton growth and 
relevance for the spring bloom. J Plankton Res 30: 75—9 1 . 

1 2 . Rhec GY, Gotham IJ (1981) The effect of environmental factors on 
phytoplankton growth: Light and the interactions of light with nitrate limitation. 
Limnol Oceanogr 26: 649-659. 

1 3 . Rhee GY, Gotham IJ (1981) The effect of environmental factors on 
phytoplankton growth: Temperature and the interactions of temperature with 
nutrient limitation. Limnol Oceanogr 26: 635-648. 

14. Ahlgren Q (1988) Phosphorus as growth-regulating factor relative to other 
environmental factors in cultured algae. Hydrobiologia 170: 191-210. 

15. Cembella AF^, Antia NJ, Harrison PJ (1984) 'Fhe utilization of inorganic and 
organic phosphorus compounds as nutrients by eukar\'otic mieroalgac: A 
multidisciplinary perspective: Part 1. Crc Cr Rev Microbiol 10: 317-391. 

16. Goldman JC (1977) Temperature effects on phytoplankton growth in continuous 
culture. Limnol Oceanogr 22: 932-936. 

17. Riegman R, Stolte W, Noordeloos AAM, Slezak D (2000) Nutrient uptake and 
alkaline phosphatase (EC 3: 1 :3: 1) activity of Emiliania huxleyi (Prymnesiophyceae) 
during growth under N and P limitation in continuous cultures. J Phycol 36: 87— 
96. 

18. Litchman E, Klausmeicr CA, Bossard P (2004) Phytoplankton nutrient 
competition under dynamic light regimes. Limnol Oceanogr 49: 1457—1462. 

19. Dickman EM, Vanni MJ, Horgan MJ (2006) Interactive effects of light and 
nutrients on phytoplankton stoichiometry. Oecologia 149: 676-689. 



adaptation. The influence of the photoperiod seems to be 
restricted to nutrient replete rather than P-limited growth rates, 
although the same need not apply to N-limitation. The nature of 

the temperature interactions with P-quota suggests that warming 
should counteract de-eutrophication and, when phosphorus is 
limiting, favour warm-adapted species more than expected based 
only on their temperature optima. 

Acknowledgments 

We thank Marion KaulfiiB for help in algae cultivation, Antje Liider for 
phosphorus measurements, Marianne Graupe for carbon measurements, 
as well as Dr. Katrin Teubner and two other anon\mous reviewers for 
helpful comments which improved the manuscript. 

Author Contributions 

Conceived and designed the experiments: AN TS JK. Performed the 
experiments: AN TS. Analyzed the data: TS AN. Contributed reagents/ 
materials/analysis tools: JK AN. Wrote the paper: TS. 



20. Litchman E, Steiner D, Bossard P (2003) Photosynthetic and growth responses of 
three freshwater algae to phosphorus limitation and daylength. Freshwater Biol 
48: 2141-2148. 

21. Riegman R, Mur LR (1985) Effects of photoperiodicity and light irradiance on 
phosphate-limited Oscillatoria agardhti in ehemostat cultures. 2. Phosphate uptake 
and growth. Arch Microbiol 142: 12— To. 

22. Castenholz RW (1964) The effect of daylength and light intensitv on the growth 
of littoral marine diatoms in culture. Physiol Plant 17: 951-963. 

23. Paasche E (1968) Marine ]:)lankton algae grown with light-dark cycles. 2. Ditylum 
brighiwellii and Mtzschia turgidula. Physiol Plant 21: 66-77. 

24. Foy RH, Gibson CE, Smidi RV (1976) The influence of daylength, %ht 
intensity and temperature on the growth rates of planktonic blue-green algae. Br 
PhycolJ 15: 151-163. 

25. Gibson CE, Foy RH (1983) The photosynthesis and growth efficiency of a 
planktonic blue-green a!ga, Osfdlaloria redekei. Br Ph\ToIJ 18: 39-45. 

26. Droop MR (1968; X'ilainin XW'I and marine ecology. 1. Kinetics of uptake, 
growth and inhibition in Alonochijsis lullieri. } Mar Biol Assoc UK 48: 689—733. 

27. Flynn KJ (2008) Use, abuse, misconceptions and insights from quota models - 
The Droop cell quota model 40 years on. Oceanogr Mar Biol 46: 1-23. 

28. Klausmeier CA, Litchman E, Daufresne T, Levin SA (2008) Phytoplankton 
stoichiometry. Ecol Res 23: 479-485. 

29. Azad HS, Borchard JA (1970) Variations in phosphorus uptake by algae. 
Environ Sci Technol 4: 737-743. 

30. Goldman JC (1979) Temperature eflects on steady-state growth, phosphorus 
uptake, and the chemical composition of a marine phytoplankter. Microb Ecol 
5: 153-166. 

31. Cembella AD, Antia NJ, Harrison PJ (1984) 'Fhe utilization of inorganic and 
organic phosphorus compounds as nutrients by eukar\'otic mieroalgae: A 
multidisciplinary perspective: Part 2. Crc Cr Rev Microbiol 11: 13—81. 

32. Ahlgren G (1987) 'Femperature functions in biologN' and their application to 
algal growth constants. Oikos 49: 177-190. 

33. Wernicke P, NickJisch A (1986) Light/dark cycle and temperature - their impact 
on phosphate-limited growth of Oscillatoria redekei VAN GOOR in semicontin- 
uous culture. Int Rev GesEunten Hydrobiol 71: 297-313. 

34. van Donk E, Kilham SS (1990) Temperature effects on silicon- and phosphorus- 
Hmited growth and competitive interactions among three diatoms. J Phycol 26: 
40-50. 

35. Flynn KJ (2008) The importance of the form of the quota curve and control of 
non-limiting nutrient transport in phytoplankton models. J Plankton Res 30: 
423-438. 

36. KohlJCi, (iiersdorf K (1991) Competition abilitv of two planktic diatoms under 
different vertical light gradients, mbdng-depth and -frequencies: An experimen- 
tal approach. Vcrhandlungcn der Internationale n Vercinigung fiir Theoretische 
und Angewandte Limnologie 24: 2652-2656. 

37. Nieklisch A (1998) Growth and light absorption of some planktonic 
cyanobacteria, diatoms and Chlorophyceae under simulated natural l^ht 
fluctuations. J Plankton Res 20: 105-119. 

38. Shatwell T, Nieklisch A, Kohler J (2012) Temperature and photoperiod effects 
on phytoplankton growing under simulated mixed layer light fluctuations. 
Limnol Oceanogr 57: 541-553. 

39. Shatwell T, Kohler J, Nieklisch A (2013) 'Femperature and photoperiod 
interactions with silicon-limited growth and competition of two diatoms. 
J Plankton Res 35: 957-971. 



PLOS ONE I www.plosone.org 



14 



July 2014 I Volume 9 | Issue 7 | e102367 



Factor Interactions in Diatoms under P-Limitation 



40. Rololi'B, Nicklisch A (1993) Partitioning of phosphate between bhie-green algae 
and their accompanying bacteria in phosphate -hmited culture. Arch Hydrobiol 
126: 405^16. 

41. Falkner G, Wagner F, SmaU JV, Falkner R (1995) Influence of fluctuating 
phosphate supply on the regulation of phosphate uptake by the blue-green alga 
Anacystis nidulans. J Phycol 31: 745-753. 

42. Healey FP (1980) Slope of the Monod equation as an indicator of advantage in 
nutrient competition. Mierob Ecol 5: 281 286. 

43. Nicklisch A (1999) Competition between the cyaiiobactcrium Limnothrix redekei 
and some s]:)ring s]:)eeies of diatoms under P-limitation. International Review of 
Hydrobiology 84: 233-241. 

44. GuiUard RRL, Lorenzen CJ (1972) YeUow-green algae with chlorophyllide c. 
J Phycol 8: 10-14. 

45. Gibson CE (1987) Adaptations in Oscillatoria redekei at very slow growth rates - 
changes in growth efficiency and phycobilin complement. Br Phycol J 22: 187— 

19l/ 

46. DEV (2007) Deutsche Einheitsverfahren (DEV) zur Wasser-, Abwasser- und 
Schlammuntcrsuchung. Analytical standard DIN 38405, part 21. Weinheim: 
VCH Verlagsges. mbH, Beuth Verlag (imbH. 

47. NickHsch A, Steinberg CEW (2009) RNA/protein and RNA/DNA ratios 
determined by ffow cytometry and their relationship to growth limitation of 
selected planktonic algae in culture. Eur J Phycol 44: 297-308. 

48. Schreiber U, Bilger W (1993) Progress in chlorophyU fluorescence research: 
Major developments during the past years in retrospect. Progress in Botany 54: 
151-173. 

49. Euhs GW (1969) Phosphorus content and rate of growth in the diatoms Cyclotella 
nana and Thalassiosira Jluviatilis. J Phycol 5: 312—321. 

50. Flynn KJ (2002) How critical is die critical N: P ratio? J Phycol 38: 961-970. 

51. CaperonJ, Meyer J (1972) Nitrogen-limited growth of marine phytoplankton. 1. 
Changes in population characteristics with steady-state growth rate. Deep-Sea 
Research 19: 601-618. 

52. Bates DM, Watts DG (1988) Nonhnear Regression Analysis and Its Applications. 
New York: Wiley. 365 p. 

53. Burnham KP, Anderson DR (2002) Model selection and multimodel 
inference: A practical information-theoretic approach. New York: Springer- 
Verlag. 488 p. 

54. R Core Team (2013) R: A language and environment for statistical computing. 
Vienna: R Foundation for Statistical Computing. 

55. Soetaert K, Petzoldt T, Setzer RW (2010) Solving Differential Equations in R: 
Package deSoIve. J Stat Softw 33 (9): 1-25. 

56. Soetaert K, Petzoldt T (2010) Inverse Modelhng, Sensitivity and Monte Carlo 
Analysis in R Using Package EME. J Stat Softw 33 (3): 1-28. 

57. Button DK (1978) On the theory of control of microbial growth kinetics by 
limiting nutrient concentrations. Deep-Sea Research 25: 1 163-1 177. 

58. LehmanJT, Botkin DB, Likens GE (1975) The assumptions and rationales of a 
computer model of phytoplankton population dynamics. Linmol Oceanogr 20: 
343-364. 



59. Tilman D (1982) Resource Competition and Community Structure. Princeton, 
NJ.: Princeton University Press. 296 p. 

60. Brown EJ, Button DK (1979) Phosphate-limited growth kinetics of Seknastnan 
capricomutiim. ] V'hycol 15: 305-311. 

61. Sterner RW (1995) Elemental stoichiometry of species in ecosystems. In: Jones 
CG, Lawton JH, editors. Linking Species and Ecosystems. New York: Chapman 
& HaU. 240-252. 

62. Darley WM (1977) Biochemical composition. In: Werner D, editor. The Biology 

of Diatoms. London: Blackwell Scientific Publications. 198-223. 

63. Sterner RW, ElserJJ (2002) Ecological Stoichiometry: The Biology of Elements 
from Molecules to the Biosphere. Princeton: Princeton University Press. 584 p. 

64. Tempest DW, flunterJR (1965) The influence of temperature and pH value on 
the maeromolecular composition of magnesium-limited and glyeerol-limited 
Aerohacter aerogene.s growing in a ehemostat. J (ien Microbiol 41: 267-273. 

65. Klausmeicr CA, Litchman E, Levin SA (2004) Phytoplankton growth and 
stoichiometry under multiple nutrient limitation. Linmol Oceanogr 49: 1463— 
1470. 

66. Giersdorf K (1988) Der Einflup von Licht und Temperatur auf das Wachstum 
der planktischen Diatomeen Stephmodiscus minutulus (Kutz.) Cleve & MsUer und 

.Nitzschia acicularis W. Smith. Mathematisch-Naturwissenschaftliche Fakultat. 
Berlin: Humboldt-Universitat zu Berlin. 117. 

67. Healey EP (1973) Characteristics of phosphorus deficiency in Anahaena. J Phycol 
9: 383-394. 

68. Lynn SG, Kilham SS, Kreeger DA, Interlandi SJ (2000) Effect of nutrient 
avaflability on the biochemical and elemental stoichiometry in the freshwater 
diatom Stephanodiscus minutulus (Bacillariophyceae). J Phycol 36: 510—522. 

69. Olsen Y (1989) Evaluation of competitive ability of Staurastrum leutkemuellerii 
(Chlorophyceae) and Microcystis aeruginosa (Cyanophyccac) under P limitation. 
J Phycol 25: 486-499. 

70. Andersen T (1997) Pelagic Nutrient Cycles: Herbivores as Sources and Sinks. 
Berlin: Springer- Verlag. 280 p. 

71. Gotham IJ, Rhee GY (1981) Comparative kinetic studies of phosphate-limited 
growth and phosphate uptake in phytoplankton in continuous culture. J Phycol 
17: 257-265. 

72. Perry MJ (1976) Phosphate utilization by an oceanic diatom in phosphorus- 
limited ehemostat culture and in oligotrophic waters of the central North Pacific. 
Limnol Oceanogr 21: 88-107. 

73. Grover JP (1989) Phosphorus-dependent growth kinetics of 1 1 species of 
freshwater algae. Limnol Oceanogr 34: 341—348. 

74. Dueobu 11, lliiisman J, Jonker RR, Mur LR (1998) (^ompelilion between a 
prochiorophyte and a cyanobacterium under various phosphorus regimes: 
Comparison with the Droop model. J Phycol 34: 467^76. 

75. Sommer U (1994) The impact of light intensity and daylength on silicate and 
nitrate competition among marine phytoplankton. Limnol Oceanogr 39: 1680- 
1688. 



PLCS ONE I www.plosone.org 



15 



July 2014 I Volume 9 | Issue 7 | e102367 



